Spontaneous generation of chromomagnetic 
fields at finite temperature in the SU(3) 
gluodynamics on a lattice 

O ■ Vadim Demchik* Alexey Gulov^and Natalia Kolomoyets^ 

^ Dnipropetrovsk National University, 49010 Dnipropetrovsk, Ukraine 

Q : December 27, 2012 

(N 

. Abstract 

The spontaneous generation of homogeneous chromomagnetic 
O ■ fields in the lattice SU(3) gluodynamics is investigated in the decon- 

finement phase of the model. A new approach based on direct mea- 
surements of the field strength on a lattice is developed. Vacuum 
magnetization is established by its infiuence on the probability den- 
XjT^ • sity function of the simulated field strength. It is found that both the 

OO I chromomagnetic fields corresponding to the diagonal SU(3) generators 

are simultaneously condensated and appear to be spatially co-directed. 
No vacuum magnetization is detected for the other SU(3) components. 
The temperature dependence of the spontaneously generated fields in 
CN \ physical units is fitted in the temperature interval 200 MeV - 200 GeV 

as the usual power law with the anomalous dimension. 

> : 
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H ■ 1 Introduction 
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The origin of cosmic magnetic fields in the early universe is one of the most 
actual problems of high energy physics. Numerous theoretical descriptions of 
large-scale magnetic field generation are discussed in the literature [HE]- One 
of the elegant mechanisms is the spontaneous vacuum magnetization at high 
temperature. The possibility of the phenomenon has already been shown 
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by both the analytical and numerical methods in the SU(2) gluodynamics 
[SI in El E], in the SU(3) gluodynamics [7], in the standard model (SM) [S], and 
the minimal supersymmetric standard model (MSSM) [9]. The stability of 
the spontaneously generated fields at high temperature was shown in Ref. |5]. 
Actually, the spontaneous vacuum magnetization is one of the distinguishable 
features of asymptotically free theories [SI [10] . 

In the SU(2) gluodynamics only one neutral magnetic field can be sponta- 
neously generated, whereas in the SU(3) gluodynamics there are two possible 
spontaneously generated neutral chromomagnetic fields [H^ and in case 
of the standard Gell-Mann matrices). The spontaneous vacuum magnetiza- 
tion phenomenon was investigated for the SU(2) gluodynamics on a lattice 
in [0], where the magnetic field was introduced with the twisted boundary 
conditions. The magnetized state was obtained by comparison of the free 
energies for the trivial vacuum state {He = 0) and the magnetized vacuum 
state {He 7^ 0). In this paper we propose a new approach to study the effect. 
In contrast to setting external fields on a lattice, it should be reasonable to 
try to detect the presence of the condensed fields directly in the standard 
Monte Carlo (MC) techniques. 

In [7] the spontaneous generation of neutral SU(3) fields in the high- 
temperature limit T ^ H^ was studied by analytical perturbative methods 
and the possibility of vacuum magnetization was shown. Nevertheless, the 
question about simultaneous generation of neutral chromomagnetic fields in 
the SU(3) gluodynamics still remains. In particular, the mutual spatial ori- 
entation of these fields should be qualified. 

In the present paper the spontaneous vacuum magnetization at high tem- 
perature is investigated in the SU(3) gluodynamics on a lattice. We devel- 
oped the procedure for detection of the condensate fields by direct analysis of 
lattice gauge configurations. This procedure is general for the deconfinement 
phase of any SU(A^) lattice gauge theory. It is grounded on the statistical 
analysis of chromomagnetic fields extracted from lattice configurations. The 
Cartesian components of quantum chromomagnetic fields have to be Maxwell 
distributed in simulations, whereas the presence of condensate fields should 
change this distribution. Relative spatial direction of the condensed chromo- 
magnetic fields is a priori unknown and can be treated as a random quantity. 
The complete electromagnetic tensor can be directly measured by taking the 
trace of a plaquette multiplied additionally by the gauge group generators. 
The deconfinement phase is spatially homogeneous, so the field tensor can 
be averaged over a lattice configuration to estimate the homogeneous vac- 
uum state beyond quantum fiuctuations. The statistical distribution of these 
averages is comparing with the Maxwell distribution to detect possible con- 
densates. We show that the vacuum is magnetized for both the neutral fields 
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if 3 s-iid Hg whereas for the other fields the vacuum magnetization is absent. 

The paper is organized as follows. Sect. 2 provides the theoretical back- 
ground of our approach to detect the condensed chromomagnetic fields. The 
results of Monte-Carlo simulations are described in Sect. 3. The temperature 
behavior of the condensate fields is derived and compared with known results 
from the literature in Sect. 4. The discussion and possible applications are 
provided in Sect. 5. 

2 Detection of the chromomagnetic fields on 
a lattice 

We consider SU(A^) gauge theory on a lattice with the standard Wilson action 
^^ = ^EE^"Tr(l-f/,.(x)), (1) 

x£A fj,<u 

where /5 = 2N / g"^ is the inverse coupling, 1 denotes the N x N identity 
matrix, and U^v{x) is the plaquette variable. The sums run over all the lattice 
sites X and directions /x, z/, and the isotropic lattice spacing a is supposed. 
The plaquette is constructed from the lattice field variables - links, 

f/^(x)=exp(mA''4(x)), (2) 

where is the 6-th SU(iV) group generator, and Ax) is the field potential 
in the continuous limit. 

Links ([2]) are subject to MC simulations. We implement standard bound- 
ary conditions: the sites on the opposite lattice edges are simply looped by 
the proper links. The result of a MC run is a Boltzmann sequence of config- 
urations with some numeric values for links. Then, any physical observable 
can be measured for each configuration and averaged over the Boltzmann 
ensemble. 

We are interested to detect a kind of Savvidy vacuum in the deconfine- 
ment phase of the SU(3) gluodynamics. This means possible existence of 
a homogeneous field strength in the whole volume being a background of 
quantum fluctuations. In contrast to the conflnement phase, which can be 
spatially fragmented to bound states, the deconfinement phase always oc- 
cupy the whole lattice, so the homogeneous vacuum can be extracted as the 
field strength averaged over all the lattice sites. Actually, the huge lattices 
can show some domain structure, but relatively small lattices are expected 
to produce the condensate in some random direction. 
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The first step to detect a non-trivial vacuum state is the measurement of 
the field strength by a given lattice configuration. In fact, the field strength 
enters the plaquette variable used to construct the Wilson action ([1]). Let us 
recall the plaquette definition: 

U^u{x) = U^{x)Uy{x + afi)Ul{x + av)Ul{x), (3) 

where fi is the unit vector in direction /i. The expansion of link variables ([2]) 
for small lattice spacing is well known (see for example [TTj): 

U,^{x) = 1 + ta'F%{x)\' - \a'Fl^{x)Fl{x)\'\'' + 0{a'). (4) 

By applying the trace, one immediately obtains the field strength squared, 
which gives the correct continuous limit for the Wilson action. 

As far as the Wilson action is considered as a good approach to field 
strength squared, we can use the same plaquette to extract the field strength 
(its flux through the plaquette, actually). This could be done by multiplying 
plaquette dlj by generators A^, then taking the trace: 



Fl,{x) = -zTt [U,,{x)X'] + 0{a'). (5) 



We will concentrate on the chromomagnetic part of the field tensor Fy^ = , 
F = —H F = H 

The average of the chromomagnetic components of ([5]) over a simulated 
lattice configuration gives us a direct measurement of the homogeneous chro- 
momagnetic field generated spontaneously in the deconfinement phase, which 
will be called the condensate for brevity. In what follows we will mark these 
averaged values as H, where vector means spatial direction. These values 
are distributed within the Boltzmann ensemble with some probability den- 
sity function p{H) (p.d.f.). We suppose Cartesian components of H to be 
Gaussian with the variance o"^, which is usual in MC simulations. The vari- 
ance can be easily measured as the sample variance taken by the MC run. 
Then, the simulated field strength H is described by the following p.d.f., 

where He is the condensate. The spatial direction of He appears to be un- 
known at each configuration generated in the MC run. It is random (at least, 
it has a random part), so p.d.f. ([6]) could be written accurately as p{H, He). 
This is the reason against naive average of the Cartesian components of H 
over the run as a way to extract the condensate value from the generated 
Boltzmann ensemble. 
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Actually we are not interested in the random direction of i^c? the question 
is what the absolute value of the condensate. Using the spherical coordinates, 
we can rewrite probability in terms of the absolute values of the vectors, H 
and He, as 



p{H)dH,dHydH, = ---L_e-(^'+^'-2^=^™^^)/2'^'i72^i7#dcos0, (7) 

where 6 is the angle between H and and is the azimuthal angle. Let 
us notice that the absolute value of the condensate He is not random. 

It is a well known fact, that the cosine of the angle between two random 
spatial vectors is distributed uniformly. It allows us to integrate ([7]) over 
angular variables obtaining the p.d.f. for the absolute value of the field 
strength: 

As a result, we get the value of the condensate He as a parameter of the 
distribution of some directly measured variable. Thus, we can determine He 
by analysis of the MC statistics obtained for H. 

In case of condensate absence. He = and p.d.f. ([8]) becomes the ordinary 
Maxwell distribution. A non-Maxwell shape of the p.d.f. means condensate 
generation. Since the p.d.f. contains only one unknown parameter He, it can 
be determined by fitting the theoretical expectation of the mean value by the 
corresponding sample mean from the MC run. 

For qualitative analysis it is convenient to rewrite Eq. (j8]) in dimensionless 
(normalized) terms. Let us introduce the quantity 

^0 = (9) 



which is just the mean value of H in case of He = 0. Then, we can consider 
dimensionless variable t] with the following p.d.f.: 

V=^, p(^) = ;^e-CVVVAsinh^, (10) 

where C = He-\/2/a depends on the condensate. 
The mean value of rj, 

oo 

16 g-cV4 f ^ V^'^V- sinh ^ rfr/ = /(C), (11) 
J 





C7r3/2 
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Figure 1: Function /(C) from Eq. (jlip used to fit the p.d.f. of the field strength. 
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Figure 2: The probability density function of the normalized field strength (Hs/Hq) 
obtained on a 18x24^ lattice for _ffg field at /3 = 8. The histogram corresponds to the lattice 
data (10 000 measurements), the dashed line describes the absence of the condensate, and 
the solid line represents the p.d.f. fit assuming the condensate existence. The field 
behaves similarly to this plot. 
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Figure 3: The probability density function of the normalized field strength (Hi/Hq) 
obtained on a 18 x 24^ lattice for Hi field at /3 = 8. The histogram corresponds to the 
lattice data (10 000 measurements), the dashed line describes the theoretical expectation 
for the absence of the condensate, and the solid line represents the fitted p.d.f. which 
confirms the condensate absence. Fields H2, H4, H^, Hq and behave similarly - they 
demonstrate the absence of condensates. 

can be associated with its estimator from lattice simulations. To compute 
the condensate, f{() must be inverted numerically. Function f{() is plotted 
in Fig. 1. 

In Figs. 2 and 3 we show an example of the described approach. The 
actual p.d.f. of the field strength is compared with the Maxwell distribution 
corresponding to the absence of the condensate {He = or ^ = 0, dashed 
lines). The parameter ( is fitted by the sample mean obtaining either zero or 
non-zero value. We also check that the fitted value of ( produces a good fit 
of the actual p.d.f. in case of condensate existence (solid lines correspond to 
(fTOj) ). As it can be seen, the neutral SU(3) component has the condensate 
whereas there is no condensate for the charged component Hi. It should be 
noted that the condensation is also present for the neutral SU(3) field H3 
and absent for all the charged fields H2, H4, H^, Hq and H^, as it can be 
expected from the theoretical reasons. 

Thus, to find the spontaneously generated field strength one should mea- 
sure the field tensor components F^y, Fy^ and F^z in lattice simulations, 
average them over configurations [H^, Hy, H^), calculate their variance over 
the run (cr^), compute the absolute values {H = ^H'^ + Hy + iJ|) and their 
mean value over the run {j] = H/Hq). Then, solving Eq. f lTTjl numerically 
with respect to C, the condensate value H^. = a(^/\/2 can be derived. 
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It should be noted that the developed method is applicable for SU(A^) 
gauge group for any in the deconfinement phase. 

3 Monte- Carlo simulations 

The main aim of the paper is the investigation of chromomagnetic field con- 
densates in the SU(3) gluodynamics. There are two neutral chromomagnetic 
fields in the model, H3 and Hg, which correspond to the diagonal Gell-Mann 
matrices. Their condensates are expected to be generated at high tempera- 
tures. 

The simulations are performed on the lattices of size up to 6 x 32^ and 
18 X 24^ for various values of /3 corresponding to the deconfinement regime. 
The lattice is updated by the heat-bath algorithm. To thermalize the system, 
1000 MC sweeps are used. All the measurements are performed on 10000 
configurations separated by 10 sweeps. The MC kernel tries to update every 
lattice link on each sweep up to 10 times. 

MC simulations are performed on graphics processing units (GPU) on 
a HGPU cluster (AMD Radeon HD 7970, HD 6970, HD 5870 and NVidia 
GeForce GTX 560 Ti). The trivial parallelization scheme is implemented - 
each GPU device simulates one lattice for a single parameter set per run. 
The subroutines for MC simulations on GPU (GPU kernels) are written in 
the Open Computing Language (OpenCL) [12]. AH the MC simulations and 
measurements are performed with the full double precision. All the necessary 
data are stored on GPU memory during a run, so the lattice size is restricted 
by the GPU device memory. 

The Cartesian components of the chromomagnetic part of the field 
strength are measured at every working configuration for all the SU(3)- 
directions. Their averages over the lattice (denoted as Ha^i, a = 1, . . . , 8; 
i = x,y, z) play a role of primary data for further statistical analysis. The 
primary data are processed for every component separately in accordance 
with the approach described in Sect. 2. First, we compute the variance cxj over 
the run for each Cartesian component Hi. By averaging these variances over 
spatial directions we obtain the variance entering p.d.f. (jS]), o" = 1/3^0",. 
Then, the value of Hq from Eq. ([9]) can be found, which is used as a unit 
for H. Finally, we calculate the normalized absolute values r] = H/Hq = 
HI + if 2 _|_ 1 ff^ appeared in the run and analyze the distribution of the 
obtained sample. 

To fit distribution (fTUj) describing the condensate generation, we com- 
pute the average over the run (sample) r/, and find the parameter ( as the 
numerical solution of Eq. (ITT]) . After that, the histogram of the sample of 77 



8 



is plotted to be compared with the fitted distribution flTUl) and the Maxwell 
distribution (the case of C = meaning the condensate absence) . An example 
of the histogram is shown in Fig. 1 for the lattice 18 x 24^ and /3 = 8 for the 
field Hs- The dashed line in Fig. 1 corresponds to the Maxwell distribution 
(no condensate), and the solid line represents the fitted distribution with the 
condensate. It can be seen that the histogram indicates the presence of the 
condensate rather then the condensate absence. Thus, the fitted parameter 
C gives the estimate of the field condensate. He = (j(/\/2. 

Analyzing all the SU(3) components of the chromomagnetic field tensor 
we find the condensates of neutral fields {H^ and Hg), On the other hand, 
no condensates are detected for the rest of the SU(3) components. 

In order to convert the obtained lattice values into the physical ones the 
standard relations are used. 



H=^, T=^, (12) 
aLf 



where ip is the flux of the condensate field through the plaquette, a is the 
lattice spacing and Lt is the lattice size in the temporal direction. The values 
of the lattice spacing a{g) have been computed using the standard technique 
described, for example, in [13]. 

The numerical results of Monte Carlo simulations for lattices 6 x 32'^, 
12 X 24^ and 18 x 24^^ are presented in Table 1. In the first column of the 
table the inverse coupling f3 is shown. The lattice spacing a in physical units 
is given in the second column. The third column contains the temperatures 
T in GeVs. The fluxes of condensate chromomagnetic fields through the 
plaquette are shown in the fourth and fifth columns. When all the quantities 
are measured in the energy units, the fluxes are dimensionless. The last 
two columns contain corresponding values of the condensate chromomagnetic 
fields if 3c and Hgc in GeV^. 

Another interesting question is the mutual spatial orientation of the con- 
densate fields and ifg. The relative angle between the third and eighth 
fields can be found for every measurement in the run, since we know all the 
Cartesian components of the corresponding vectors. The cosine of the angle 
is 

cos^ = Tfrr^ • (13) 

In Fig. 4 we show the histogram for cos 6'* obtained for lattice 18 x 24^ at 
(3 = 8. The fields seem to be co-directed, since the p.d.f. does not evi- 
dently fit the uniform distribution expected for two non-correlated vectors. 
To provide statistical estimates, let us note that random quantity f|T3l) is ex- 
pected to be distributed as the sample correlation for two Gaussian quanti- 
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Table 1: The condensates of chromomagnetic fields for different lattices. 



/3 


a, 














xlO-2 fm 


GeV 


xlO-3 


xlO-3 


GeV^ 


GeV^ 


12 X 


8 


1.994 


0.825 


3.716 


2.239 


0.7277 


0.4385 


9 


0.6335 


2.596 


3.056 


1.781 


5.931 


3.456 


10 


0.2002 


8.213 


2.582 


1.487 


50.17 


28.89 


11 


0.06301 


26.10 


2.282 


1.325 


447.7 


259.9 


12 


0.01976 


83.22 


1.931 


1.160 


3853 


2314 


18 X 243 


8 


1.994 


0.550 


3.710 


2.176 


0.7266 


0.4261 


9 


0.6335 


1.730 


3.108 


1.797 


6.032 


3.486 


10 


0.2002 


5.476 


2.649 


1.589 


51.47 


30.88 


11 


0.06301 


17.40 


2.241 


1.299 


439.7 


254.9 


12 


0.01976 


55.48 


1.987 


1.144 


3963 


2283 


6 X 32^ 


8 


1.994 


1.649 


3.686 


2.239 


0.7219 


0.4384 


9 


0.6335 


5.191 


3.113 


1.809 


6.040 


3.511 


10 


0.2002 


16.43 


2.610 


1.545 


50.71 


30.02 


11 


0.06301 


52.20 


2.190 


1.302 


429.7 


255.3 


12 


0.01976 


166.4 


1.863 


1.105 


3716 


2204 



ties. Namely, applying Fisher's transformation, we obtain Gaussian quantity 
artanli(cos^*), which mean and variance can be easily computed. In Fig. 4 
we also present the fitted distribution, the mean value of artanh(cos ^*) cor- 
responds to cos 9* ~ 0.829 confirming the hypothesis of co-directed fields. 

4 Temperature dependence of the field con- 
densates 

In the previous section it was found that both the chromomagnetic fields H^c 
and Hsc are generated in the SU(3) gluodynamics (see Table 1). It can be 
seen that the field condensate values are temperature dependent. In order to 
analyze this dependence, we plot the field strength versus the corresponding 
temperature squared. These plots for H^c and i/gc for lattices 6 x 32"^, 8 x 28'^, 
12 X 24^ and 18 x 24^ are shown in Fig. 5 in the logarithmic scale. The field 
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Figure 4: Probability density function of the cosine of the angle between the condensates 
H^c and Hsc on 18 x 24'^ lattice at /3 — 8. The solid line represents the fit using Fisher's 
transformation. 

strengths for various temperatures for each lattice size are plotted by different 
markers. The scatter plots demonstrate linear dependence in logarithmic 
coordinates for various lattice sizes. The field condensates are practically 
independent of the lattice size in spatial and temporal directions for fixed T, 
so finite-size effects seem to play no crucial role in the observed effect. 

Table 2: Temperature dependence of the condensate fields fitted by the power law H = 
fcT2(T/To)-^ with To = 0.2 GeV. 









lattice 


k 


b 


k 


b 


6 X 32^ 


0.33 


0.12 


0.20 


0.12 


8 X 28^ 


0.60 


0.13 


0.38 


0.14 


12 X 24^ 


1.29 


0.13 


0.75 


0.13 


18 X 24^ 


2.65 


0.13 


1.52 


0.12 



Obviously, the plotted temperature dependence of the field condensates 
can be fitted by the power law 

H, = kT'(^^^ \ (14) 
where Tq is some arbitrary normalization scale. We choose Tq = 0.2 GeV. 
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Figure 5: T^-dependence of the field strength H-^c (left) and iJsc (right). The dashed 
lines correspond to the data fits by Eq. (fT4)) . The solid line corresponds to the analytic 
perturbative approach [7]. 



The fit results for eacli lattice size are sliown in Table 2. Tlie first column 
contains the information about lattices, the same as in Fig. 5. The fitted 
parameters k and h for the field are present in the second and third 
columns. The same parameters corresponding to the field are shown in 
the last two columns. As it can be seen from Table 2, the slope parameters 
h are very close for different lattices. This parameter is just 'the anomalous 
dimension' of the field condensate, and it can be set to 0.13 for the SU(3) 
gluodynamics. The coefficient k, which defines the level of the fitting curve, 
grows with the increasing of temporal lattice size, but remains independent 
of the spatial one. This means, probably, either some finite-size effects, or 
systematic errors of the approach. 

As it can be seen from above, both of the fields has similar magnitudes and 
identically grow in the investigated temperature range. No field dominates 
in the temperature range considered. 

As it can be expected from the dimensional analysis, the dependence of 
the condensate field strength on the temperature is close to the quadratic 
function, if ~ T^. The non-trivial behavior is described by the phenomeno- 
logical law with the anomalous dimension h. It is interesting to compare the 
obtained results with analytic approaches assuming the continuous space- 
time. In [7] the generation of chromomagnetic fields at high-temperature 
limit (T 3> y/gH) was derived within the perturbative effective potential 
including the one-loop and daisy diagrams. Therein the following law was 
claimed 

Hs, = 0.2976^, Hsc = 1.8351^, (15) 
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where g = g{T) is the running couphng. We can take the QCD running 
couphng to estimate the temperature behavior of H. At high temperatures 
the QCD couphng constant a = g"^ / (47r) can be written as 



1.38808 



(16) 



a = 



ln(T2/(0.217 GeV)2)' 



where the numeric parameters were obtained by fitting the data from [H]. 
The temperature dependence of H3 and Hg from f|T5|) is shown in Fig. 5 by 
the sohd hues. The slope and level of these lines are in a good agreement 
with the phenomenological law f[T^ from our MC results. 

5 Discussion 

In the paper spontaneous generation of chromomagnetic fields in the lattice 
SU(3) gluodynamics is investigated in the deconfinement phase. An approach 
for the detection of homogeneous condensate fields by direct measurements 
of the Cartesian components of the field tensor on a lattice is developed. The 
method is based on two steps. First, the homogeneous background field is 
found at each simulated configuration by averaging the components of the 
field tensor over the lattice. Second, the presence of condensate fields affect 
the statistical distribution of the extracted homogeneous background over the 
MC run. The p.d.f. of the absolute value of the background field depends only 
on the absolute value of the condensate field being distributed in accordance 
with the Maxwell distribution in case of no condensate fields. A non-Maxwell 
shape of the distribution indicates spontaneous generation of the condensate 
fields and allows to fit the value of the generated condensate. Probing differ- 
ent lattice geometries and various temperatures (see Table 1 and Fig. 5), we 
found the evident spontaneously generated magnetic components of the field 
tensor in the 'neutral' SU(3) directions (both in the third and eight ones) 
whereas other SU(3) directions show the clear zero background. This serves 
as an additional argument that our approach refiects well motivated physics 
and cannot be interpreted as some random artificial lattice effects. In the 
physical units, the strength of condensate chromomagnetic fields monoton- 
ically grows with temperature increasing. We successfully fit the collected 
MC data by the phenomenological law with 'anomalous dimension': 



This law incorporates the results obtained by analytic perturbative approach 
[7] as a high-temperature asymptotics. 




(17) 
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Both the neutral chromomagnetic fields H^c and ifgc are simultaneously 
generated in the deconfinement phase in the whole temperature range in- 
vestigated. These condensate fields have the same functional temperature 
dependence, the same order of magnitudes, so they evolve in a balanced 
manner. The direct analysis of the Cartesian components of the condensed 
fields shows that the fields are spontaneously generated in the same spatial 
direction. 

The main results of the paper are obtained in physical units. It is checked 
that the condensate field strengths are independent of the spatial lattice size. 

In general, the homogeneous chromomagnetic field in the whole space 
breaks the gauge invariance. Thus, the formation of domain structure at 
large scales is expected, which should restore the gauge invariance. So, the 
investigation of domain structure on huge lattices is also an interesting sub- 
ject of further studies. 

The method developed in Sect. 2 to determine the magnetized vacuum 
state is quite general. It can be applied for any SU(A^) lattice gauge theory. 
In this paper we restrict the investigation to the SU(3) gluodynamics as an 
important example. 

In Sect. 4 we show that both the chromomagnetic fields H^c and Hsc are 
simultaneously generated in the deconfinement phase in the wide tempera- 
ture interval from 200 MeV to 200 GeV. These condensate fields have the 
same functional temperature dependence, so they evolve in a balanced man- 
ner. Moreover, both the fields appear to be spontaneously generated in the 
same spatial direction. It is derived in MC simulations that the spontaneous 
vacuum magnetization takes place only for and Hg chromomagnetic fields 
{Hsc 7^ 0, Hsc 7^ 0), whereas there is no condensation of other chromomag- 
netic fields {Hic, H2c, Hie, H^c, -f^ec Hj^ = 0). 

In [8] it was shown in the 'one-loop plus daisy' perturbative approach 
that the chromomagnetic field H^c is generated in the standard model. The 
source of spontaneous magnetization effect is in the gauge sector of the model. 
The fermionic sector of the standard model slightly weakens the effect only. 
The spontaneous vacuum magnetization was also studied in the MSSM [9], 
where the fermionic sector of the model is extended with the s-particles. 
It was obtain therein that doubled fermionic sector of the MSSM does not 
eliminate the effect. As for our results, fermions should be accounted in 
further investigations. 
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